Critical Behavior of Frustrated Josephson Junction Arrays with Bond Disorder 
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The scaling behavior of the current-voltage (IV) characteristics of a two-dimensional proximity- 
coupled Josephson junction array (JJA) with quenched bond disorder was investigated for frustra- 
tions / = 1/5, 1/3, 2/5, and 1/2. For all these frustrations including 1/5 and 2/5 where a strongly 
first-order phase transition is expected in the absence of disorder, the IV characteristics exhibited 
a good scaling behavior. The critical exponent v indicates that bond disorder may drive the phase 
transitions to be continuous but not into the Ising universality class, contrary to what was observed 
in Monte Carlo simulations. The dynamic critical exponent z for JJA's was found to be only 0.60 - 
0.77. 
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The addition of bond randomness may change the first- 
order transition into the second-order one. A particular 
manifestation of the rounding of the first-order transition 
is seen pi 0] in the two-dimensional (2D) g-state Potts 
model, for which the transition in the absence of disorder 
is of first-order for q > 4. For the eight-state Potts model, 
the critical exponents of the disorder-driven second-order 
transition were found H to fall into the universality class 
of the 2D Ising model. Kardar et al. H suggested that 
the disorder-induced Ising criticality is due to the asymp- 
totically Ising-like critical interface. It has been shown in 
the recent Monte Carlo (MC) simulations f| that bond 
randomness has a similar effect on the phase transition of 
the uniformly frustrated 2D XY model. When quenched 
bond disorder was introduced, the first-order transition 
of the XY model at frustration, the fraction of a flux 
quantum per plaquette, / = 2/5 was found to change 
into the continuous transition with the critical exponents 
of the 2D Ising values. Despite the difference in the un- 
derlying symmetry of the XY model from the g-state 
Potts model, the disorder-induced Ising-like transition at 
/ = 2/5 was also attributed to the Ising-like domain wall 
excitations. The numerical observations were claimed to 
be consistent with the critical exponents obtained exper- 
imentally However, because of a systematic error 
in locating T c of the sample, a superconducting wire net- 
work, the announced error bar for the experimentally ob- 
tained critical exponents is ±0.5, somewhat too large for 
identifying the Ising universality class. It appears to us 
that an experimental confirmation of the Ising criticality 
at / = 2/5 with disorder has yet to be done. Another 
interesting question raised by the numerical observations 
at / = 2/5 is whether the transitions at different frus- 
trations such as 1/5, 1/3, and 1/2 can also be driven to 
the same Ising universality class by bond disorder. In the 
absence of disorder, expected numerically are a strongly 
first-order transition at / = 1/5 H, an Ising-like transi- 
tion at / = 1/3 0, and a continuous but non-Ising-like 
transition at / = 1/2 B. 



Proximity-coupled Josephson junction arrays (JJA's) 
and superconducting wire networks (SWN's) exposed to 
a uniform transverse magnetic field are the closest exper- 
imental realizations of the uniformly frustrated 2D XY 
model. In JJA's and SWN's, the phase transition ap- 
pears in the form of a superconducting-to-resistive tran- 
sition associated with the melting of a vortex lattice. If 
the superconducting transition is a continuous one, the 
critical behavior of the phase transition can be investi- 
gated experimentally by examining the scaling behavior 
of the current-voltage (IV) characteristics of the super- 
conducting sample |fij, |jj . For the purpose of determin- 
ing the critical exponents, JJA's gain an advantage over 
SWN's. A broad critical region of the JJA allows one 
to determine the critical exponents with accuracies suf- 
ficient for identifying the universality class. We investi- 
gated the scaling behavior of the IV characteristics of a 
square proximity-coupled JJA with the deviation of the 
junction coupling strength < 15 % at four different frus- 
trations / = 1/5, 1/3, 2/5, and 1/2. The scaled IV data 
reveal that the phase transitions at all these f's are con- 
tinuous, but none of the transitions including those at 
/ = 2/5 and 1/3 belongs to the Ising universality class. 
It was also found from the scaling analyses that the dy- 
namic critical exponent z of the JJA's is only 0.60 - 0.77, 
much smaller than that observed in SWN's. 

The sample studied consists of an array of 200x1000 
Nb/Cu/Nb Josephson junctions. The 0.2-/im-thick cross- 
shaped Nb islands patterned on top of a 0.3-/xm-thick Cu 
film form a square array with a lattice constant of 10 (im 
and the junctions are 4 //m wide and 1.3 /an long. The 
zero-field superconducting transition temperature of the 
sample T KT = 6.24 K and the width ~ 0.4 K. The de- 
viation of the single-junction critical current i c or the 
junction coupling strength J (— hi c /2e) of the sample, 
estimated from the I vs dV/dl curve in Fig. [j](a) and 
the junction separation fluctuations, was < 15 %. The 
appearance of higher-order dips in the magnetoresistance 
curve of the sample with 10-/iA excitation current, shown 
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FIG. 1: (a) / vs dV/dJ at T = 4.40 K. (b) Sample resistance 
with 10-fiA excitation current as a function of frustration at 
T= 5.13 K. 



in Fig. |l|(b), indicates the good uniformity of the ap- 
plied transverse magnetic field over the sample. From 
the magnetoresistance curve of the sample with 50-/iA 
excitation current, showing sharp resistance minima at 
fractional /'s, a precise adjustment of the magnetic field 
or the frustration was possible. The temperature during 
the measurements was controlled with fluctuations < 1 
mK. The standard four-probe technique was adopted for 
measurements of the IV characteristics. The voltage V 
was measured by the use of a transformer-coupled lock- 
in amplifier with a square-wave current at 23 Hz. The 
single-junction critical current i c at low temperatures can 
be obtained from the I vs dV/dl curves. The i c 's at high 
temperatures were determined by extrapolating the i c 
vs T data at low temperatures by the use of de Genncs 
expression for proximity-coupled junctions in the dirty 
limit, i c (T) = i c (0)(l - T/T co ) 2 exp(-aT 1 / 2 ), where T co 
is the BCS transition temperature. 

Figure || shows some of the IV traces for / = 1/5, 
1/3, 2/5, and 1/2. The IV curves at high tempera- 
tures are concave upward. As the temperature is low- 
ered, the curves become straight and then bend progres- 
sively more downward jl(| . The temperature dependence 
of the IV curves suggests that for all the /'s studied, 
the sample experiences a superconducting transition at 
the temperature where the straight IV curve appears. 
If the superconducting transition is continuous, the IV 
data are expected to satisfy the scaling relation in 2D, 
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FIG. 2: Some of the IV curves for four different frustrations, 
(a) / = 1/5 at T = 5.120 to 5.950 K; (b) / = 1/3 at T = 
5.400 to 5.999 K; (c) / = 2/5 at T = 5.249 to 5.876 K; (d) 
/ = 1/2 at T = 5.600 to 6.100 K. The dashed lines are drawn 
to show the power law (V ~ I z+1 ) behavior at the critical 
temperature. 



V/I\T - T c \ zy = £±(I/T\T - T c \"), where £± are the 
scaling functions above and below T c jnj. This scaling 
form becomes a simple power-law IV relation, V ~ I z+1 , 
at T — T c . Thus we can find approximate T c and z di- 
rectly from the IV curves in Fig. |[ Accurate T c , z, and 
v can be obtained from the scaling analyses of the IV 
data. Though the scaling analysis of IV data is a pow- 
erful tool for identifying a continuous superconducting 
transition and for determining the critical exponents at 
the transition, it has been recently shown |l2|, [l]| that 
special care should be taken not to extract any false infor- 
mation from the scaling analysis. A good collapse of IV 
curves is not sufficient evidence for the transition. The 
proposed necessary criterion for IV data supporting a su- 
perconducting transition is that log V vs log I isotherms 
equally distanced from T c must have opposite concavi- 
ties at the same applied currents |L2]. Figure || shows 
the current dependence of the derivative of the log V vs 
log I isotherms in Fig. |^. The data at high currents are 
due to the single-junction effect (/ = 1/5) or the free flux 
flow (/ = 1/3, 2/5, and 1/2) and thus should be excluded 
from scaling. The derivative plots clearly show that the 
isotherms near T c have the derivative nearly independent 
of the current and that the isotherms at temperatures 
equally distanced from T c have opposite concavities at 
the same current level. For all four /'s, our data appear 
to meet the criterion. 

The IV data scaled on the basis of the scaling form 
are shown in Fig. ||. Each plot contains TV curves at 
16-24 different temperatures. We find the four sets of 
data exhibit good scaling behaviors, which suggests the 
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FIG. 3: Current dependence of the slope of the IV curves in 
Fig. |. 



continuous phase transition for all the /'s. As discussed 
earlier, the continuous transitions at f — 1/5 and 2/5 for 
which a first-order phase transition is expected [|7|, [l4| are 
probably due to the random bond disorder in the sample. 
In the scaling analyses, the parameters T c , is, and z were 
carefully adjusted until a good collapse of the curves was 
achieved. The insets of Fig. |] exhibit the values of the pa- 
rameters used to scale the data. Table Q summarizes the 
results of the scaling analyses together with the results 
of numerical studies of the XY model without disorder 
H 0, H in parentheses for comparison. Both T c and v 
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from the IV measurements are higher than the numerical 
predictions in the absence of disorder. We note that for 
all the /'s including 2/5 and 1/3, the critical exponent 
v considerably exceeds the 2D Ising value [y = 1). It is 
interesting that the v (= 1.8) for / = 2/5 is close to the 
SWN value (1.8 - 1.9) at the near-irrational frustration 
/ = 0.618 ~ (V5- l)/2 §, §, the complementary frac- 
tion (1 — /) of which is 0.382 in the vicinity of 2/5. The 
results of critical scaling analyses of the IV data indicate 
that quenched bond disorder may drive the phase tran- 
sitions of frustrated JJA's to be continuous but not into 
the Ising universality class, contrary to the observed in 
simulations at / = 2/5 with disorder and at / = 1/3. 

It was argued |J that the Ising-like transitions of the 
XY model at / = 1/3 and at / = 2/5 with disorder are 
because of domain wall excitations with configurations 
of the Ising model. The arguments suggest that the ex- 
citation configurations of the JJA's in the IV measure- 
ments must be different from those of the XY systems 
in the simulations. The larger r/s for JJA's imply that 
the JJA has a much longer correlation length near T c 
than the corresponding XY system. For the present, we 
have no plausible explanation for the disagreement be- 
tween experiments and simulations. We speculate only 
about the disagreement. Though JJA's are known to 
provide a near-ideal implementation of the XY model, 
JJA's in the IV measurements could differ in the follow- 
ing aspects from the XY model in the simulations. In the 
experiments, the critical behavior was studied on JJA's 
with a current injected. If an injected current alters fun- 
damental ingredients of the system, the nonequilibrium 
measurements might bring out the critical exponents dif- 
ferent from those in equilibrium. On the other hand, 
in the simulations, the periodic boundary condition was 
imposed to compensate the effects of a finite size of the 
model system investigated. It has been found in MC sim- 
ulations of the 2D XY model near incommensurability 
|]l5| that when the system has a natural period different 
from the imposed period, the periodic boundary condi- 
tion induces numerous extra low-energy walls, which may 
have a significant effect on transitions in smaller systems. 
If it is also the case for / = 1/5, 1/3, 2/5, and 1/2, the 
critical exponents from the simulations at the /'s could 



TABLE I: T c , v, and z from the scaling analyses of the IV 
data. The numbers in parentheses are those from numerical 
studies of the XY model without disorder ||, [?], ^| . The values 
of T c are in unit of J/fcs. 



FIG. 4: Scaling plots of the IV curves. The values of T c , v, 
and z used to scale the data are shown in the insets. 
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be different from those of real arrays with larger sizes and 
free boundaries. 

Another remarkable feature of our data is that the dy- 
namic critical exponent z of the JJA's is only 0.60 - 0.77, 
much smaller than that (z — 2.0 ± 0.5) observed || in 
SWN's which are known in the same static universality 
class as JJA's. This implies that in the presence of a 
magnetic field, JJA's relax near T c much more rapidly 
than SWN's. The JJA values for z stand closer to the 
value from Langevin-dynamics simulations (z ps 1 for 
/ = 1/2) fT§ than to the SWN value. The Langevin- 
dynamics simulations of the resistively shunted-junction 
model and of the resistively capacitively shunted-junction 
model show that the IV characteristics of a JJA for 
/ = 1/2 are dominated by the motion of domain walls 
rather than the motion of vortices, z ~ 2 is seen at 
the transition of the frustrated XY model governed by 
the purely dissipative dynamics (without conservation of 
order parameter) |l7[] as well as at a Kosterlitz-Thouless- 
Berezinskii transition. We speculate that such low z's 
for the JJA's are possibly related to the brisk motion of 
domain walls. It is, however, not clear whether z ~ 2 for 
SWN's is from stronger dissipations or from some other 
origins. 

In summary, the scaling behavior of the IV character- 
istics of a square JJA with quenched bond disorder was 
investigated for f = 1/5, 1/3, 2/5, and 1/2. For all the 
frustrations, a good scaling behavior of the IV character- 
istics was found. The critical exponent v exceeding the 
2D Ising value suggests that bond disorder may drive the 
transitions of JJA's continuous but not to the 2D Ising 
universality class, contrary to what was observed in MC 
simulations of the XY model. The dynamic critical expo- 
nent z as low as 0.60 - 0.77 indicates that the frustrated 
JJA's relax rapidly in the vicinity of T c . 
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Ministry of Education. 
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